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Abstract 

We study the initial distribution of a parton system which is formed just after relativistic heavy 
ion collision by the elastic scattering among the constituent partons in details and analyze the 
baryon and strangeness contents of the primary parton system. We present the rapidity and 
energy distributions of the system. 
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I. INTRODUCTION 



The relativistic heavy ion collision at RHIC or LHC energy liberates lots of constituent 
partons of nuclei by primary collisions and those partons evolve according to quantum theory 
and may become a quark-gluon plasma which has been searched intensively by AGS and 
RHIC. The study of the early system of those partons thus are very important and we 
somehow have to know the momentum distribution as well as the space distribution to 
understand in details the dynamics of the system as a function of time. 

There are several schemes to obtain the initial (phase space) distribution of partons: 
The constituent partons of a nucleus can be characterized by the Feynman x where the 
soft parton has small x and the hard one large x. The soft parton can be identified as the 
non-abelian Weizsacker- Williams field radiated from the hard partons and is dominant in 
number. While the two nuclei pass through each other in collision, the high x partons of 
one nucleus can be scattered off from those high x partons of the other nucleus and the 
non-abelian Weicsacker- Williams field lose the coherence of the parent parton and become 
particles. This initial parton information has been studied by Krasnitz et al[lj and Lappi^ 
based on CGC^, 4|, namely, shattering of the color glass condensates of two colliding nuclei, 
and the analytic distribution function has been proposed. It however can be applicable to 
the central rapidity region only and the production of quarks or antiquarks can not be 
calculated from the scheme. 

Another theory to produce dynamically the partons after a heavy ion collision is that 
while the colliding nuclei is passing through each other, they rapidly decelerate and produce 
Schwinger-like radiation already in thermal spectrum^. The background field are the pulse- 
like CGC of McLerran-Venugopalan model of the other nucleus. This theory could provide 
the scenario to explain the fast thermalization of the parton system. 

Even though we can understand the gluon production from these studies, those do not 
provide the initial baryon density which is shown in RHIC experiments 6, 7]. To understand 
this phenomena we certainly have to include quarks as well as antiquarks from the primary 
collision. Consider the CM frame of two colliding nuclei at high energy. The colliding 
nuclei are Lorentz contracted to become disks and the transversal motion of constituent 
partons can be ignored since the longitudinal momentum is dominant at high energy. When 
the colliding nuclei pass through each other, the constituents of a projectile nucleus make 
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scattering with those of a target nucleus(note that the partons in the same nucleus do not 
collide with each other since the time delation) and become particles if the momentum 
transfer is large enough. This, so-called, factorization method, has been used to calculate 



jets or minijets 



0, 1, m 



from early 80s. There are many schemes to obtain the parton 



distribution of a hadron or a nucleus [13j]: for examples, we can fold the parton distribution 
of a nucleon and the nucleon distribution of a nucleus [m, The wave function of a high 
energy nucleus has been studied extensively and we use GRV98, CTEQ and an emperical 
distribution in our study. Since this method only gives relatively high p± partons but soft 
partons can not be obtained so that it gives a reasonable distribution at higher energy but 
does not give good at RHIC energy. The studies on this direction had been persuited by 



Eskola[8J], by Hammon et al. 10J, by Nayak et al.[9j] and by Cooper et al. ll|. We however 
think it is worth to extend the study to include baryon and strangeness contents at RHIC 
and LHC energy. We also briefly note here that the spatial distribution of produced partons 
can be understood assuming that the quanta becomes particle after the time 1/Et, which 
Et is the transversal energy of the quanta and at the place where the quanta travels p/EEt 
from the scattering point. 

We rewrite the basic formulas to calculate the production of partons in Section II. We 
present the numerical results in Section III. We summarize and conclude in Section IV. 

II. INITIAL PARTON PRODUCTION 

We assume that partons are produced at relativistic heavy-ion collisions by elastic scat- 
tering between the constituent of a projectile nucleus and that of a target nucleus. Since 
the scattering cross section gives the number of collision events per unit overlap area, we 



can write the total number of events 



11 



12| : 



dy,dy 4 d 2 PT *i/</a(*i. Ql)^fj/B(x2, Ql) dP"^ ( 21 ) 

ij, kl 

where a parton i of the nucleus A makes a collision with a parton j of the nucleus B and 
produces partons k and I. Each parton has rapidity yi, y 2 , y% and ?/ 4 , respectively. x\ and 
%2 are the Bjorken scaling variables of parton i and j. The relations between the variables 
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before and after the collision are given by 

Xl = PT ( e yi +e m )/^, (2.2) 

x 2 = p T {e- yi +e- y2 )/^, (2.3) 

s = Xix 2 s, (2.4) 

i = -p 2 T (l + e y2 - yi ), (2.5) 

u = -p 2 T (l + e yi - y2 ). (2.6) 

We also rewrite the available kinematic region for convenience 

Oo 2 <PT ! <(^-) 2 , (2.7) 
2 cosh y 

- log(^ - e"*) < y, < log(^ - e-v), (2.8) 
Pt Pt 



is the K-factor to include the higher-order diagrams; we will set K = 2 throughout our 
study. Qo is the momentum scale which we are looking at the nucleus and is a minimum 
momentum transfer. The hat on the Mandelstan variables means that those are the variables 
of a parton. 

T(b) is the nuclear geometric factor at a given impact parameter b. Assuming that the 
density of nucleus is constant over the sphere of radius R with the sharp edge, we can write 



T{b) = J dfT AB (f;b) (2.10) 
= 4p£p* J df^Rl - (f - b/2)^R% - (r + b/2)\ 

where Tab is the nuclear thickness function. Note that we assume there is no correlation 
between the position and the momentum of a given parton. 

Since each event gives two partons, we can write the distribution of produced partons as 
follows 

—— — = KT(b) I dy 4 —y^x l f l /A{xi,p 2 T )x 2 f j /B{x2,p 2 T)- — , (2.11) 

dp T dy J s ^ ' 2 dt 

Note also that a quark minijet of (y,pr) from a process, for example, gq — > gq, comes from 
the scattering of a quark of the nucleus A and a gluon of the nucleus B via the t-channel 
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or from the scattering of a gluon of the nucleus A and a quark of the nucleus B via the 
u-channel so that we have to include both the a gq ^ gq (t,u) and the a gq ^ gq (u,t) processes 
multiplied by an approperate distribution, respectively. 
The processes we consider in our study are 

99 *~> 99, qq, (2.12) 

gq <-»■ gq, (2.13) 

gq <-> gq, (2.14) 

q a q b <-> gV, (2.15) 

qq <-> gg, (2.16) 

g°g 6 <-> g c g d . (2.17) 

We however ignore some of basic channels such as qg — > 57, gg — > 77, gg — > 57, which could 
provide important information on the system. 

The cross sections for the processes up to the leading order (LO) are 



da 99 ^ 99 9na 2 tu su st 
"V 3 - 75. ~~ 77 ~~ 75.) 



dt 2s 2 " s 2 t 2 u 



da"^ q ^ b na 2 u i 9t 2 + u 2 , 
df~ ~ ~6P ah[ l + ~~ ~A~^' 



da 9q ^ 9q 4ira 2 u s 9 s 2 + u 2 
dt ~ ~9S 2 ~^1 ~~ 1 + 4 £2 J 



da qaq b ^qaq b AtTO 2 J 2 + V? t 2 + § 2 2 § 2 



d a q a q b ^qcq d Ana 2 S 2 + U 2 t 2 + U 2 2« 2 



s -5 ab h + ~ ~ ~ 
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and 



dt 27s 2 au 't u As 



do 9q ~* 9q da 9q ^ 9q 



dt dt 
da qq ^ qq da qq ^ qq 

dt di 



2.18) 
2.19) 
2.20) 
2.21) 
2.22) 
2.23) 
2.24) 
2.25) 
2.26) 
2.27) 
2.28) 

2.29) 
2.30) 



We consider a head-on collision {b = 0) and use the running coupling constant a s 
{ " where Q is the momentum transfer, b = 11 — |ny, and Aqqd — 200 MeV. 



The production of partons strongly depends on the minimum momentum transfer Qq. 
We will set in between Qq = 1.2 and 2.0 GeV/c at RHIC energy and Qo — 2.0 and 4.0 
GeV/c at LHC energy. These partons have quite high energy so that they will be called 
minijets in this study. Since the nucleons in a nucleus are treated independently, we can 
write the parton distribution of the nucleus A as follows; 

f l/A (x, Q 2 ) = f i/N (x, Q 2 )R A (x, Q 2 ), (2.31) 

where fi/N{x,Q 2 ) is the parton distribution within a free nucleon and Ra{x,Q 2 ) is the the 
nucleus ratio function, which is the distribution of a nucleon within the nucleus. 



We use the CTEQ4(17| or the GRV980 distribution function for a free nucleon parton 



distribution and the EKS98 parametrization for the ratio function [1SJ. We also use the 
emperical distribution of a nucleus Ji], Q], 

xG(x, Q 2 ) = A \og( Q2 ~t 2 A2 )x~ x (l - x) n (2.32) 

where we set A = 0.3, A = A = 0.2, n = 4 and A = 0.2. 
We compare our data with those of Krasnitz et al. [l|; 

dN nR 2 



dyd 2 k± g 



2 



f(k±/A s ) (2.33) 



where 



= 1_ \ a^MVW+^/Teff) - IY\ k ± /A s < 1.5 
1 9 2 | a 2 Ai\og(ATrkJA s )k^\ k ± /A s > 1.5 

and a x = 0.137, a 2 = 0.0087, m = 0.0358A,, T eff = 0.465A S 



III. RESULTS AND DISCUSSIONS 

Fig. 1 shows the total energy liberated after collision at RHIC energy. Hirano distribution 
gives less energy since the distribution does not include quarks and antiquarks. GRV and 
CTEQ function give almost the same results. The available total energy is about 39.4 TeV 
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Total Liberated Energy at RHIC 




1.2 1.4 1.6 1.8 2.0 



Q (GeV/c) 



FIG. 1: Total energy liberated from Au-Au at 200 GeV/pair. 

for the head-on collision at RHIC. There should be lots of soft partons and they will take 
substantial amount of energy. If we estimate that the total energy liberated to the minijets 



is 50-60% of total available CM energy 21], we can set the minimum momentum transfer to 
be in-between 1.8 and 2.0 GeV/c, and for 60-70% liberated energy, Q is in-between 1.7 and 
1.9 and for 70-80% about 1.7 GeV/c. These number are larger than those of theoretically 
predicted, Q$ = 1 — 2GeV 2 /c 2 . We further notice that the Hirano distribution predictes 
much lower values than those of GRV and CTEQ distribution but gives reasonable theoretical 
value for the minimum momentum transfer scale, Qq. 

Fig. 2 shows the total energy liberated at LHC which has 1140 TeV available CM energy. 
GRV simulation gives more energy deposition than the other two distributions. Assuming 
that the liberated energy into minijets is about 70-80%, we can estimate the minimum 
momentum transfer is 3.5-3.7 GeV/c for CTEQ distribution and 3.7-3.9 GeV/c for GRV 
model. These numbers are much larger than those of educated theoretical prediction, 
Ql = 2 - 3GeV 2 /c 2 . 

Fig. [3] shows the total number of produced partons at RHIC as a function of Qo and the 
number of produced partons is about 3700 with CTEQ at Qo = 1.7 GeV/c, which liberates 
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Total Liberated Energy at LHC 
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FIG. 2: Total energy liberated from Pb-Pb at 5.5 TeV/pair. 

about 70% of total energy. Fig. H] shows the total number of produced partons at LHC as 
function of Qq. The total number of partons with CTEQ is about 11700 at Q = 3.5 GeV/c 
which liberates about 80% of CM energy at LHC. CTEQ predicts more parton production 
than the others at RHIC but GRV gives far more partons at LHC, which means GRV 
distribution have more increasement than the others as x gets smaller. 

Table 1 gives the total numbers and the average energy in GeV of partons after a Au-Au 
collision at 200A GeV and a Pb-Pb collision at 5.5A TeV. The table shows that most of 
partons produced at LHC energy are jets and the liberated valence quarks have extremely 
high energy. 30% of valence quarks are liberated from the colliding nuclei at RHIC and 32% 
at LHC so that substantial amount of baryons are provided into the parton system from the 
very beginning. We also note that the strangeness is produced from the primary collisions. 
These initial baryon and stageness contents should be imprinted in further evolution of the 
system. 

Fig. [5] shows the number of gluons produced and energy distribution as a function of 
rapidity at RHIC energy. 



S 




FIG. 3: Total number of partons produced at RHIC energy. 




FIG. 4: Total number of partons produced at LHC energy. 

Even though the gluons are flat and high at central rapidity region, those gluons are strongly 
oriented in forward and backward direction energetically. We estimate the j^/^ ~ 2.3 
GeV at central rapidity region. This average energy of parton at central rapidity region is 
much lower than that of minijets. Fig. [6] shows the number of u-quarks produced and the 
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RHIC: Q = 1.7 GeV 


LHC: Q = 3.5 GeV 
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14.5 


384 


253.3 
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33 


7.9 


122 


79.9 
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62 


8.1 


170 


90.2 
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61 


8.2 


168 


92.0 



TABLE I: Initial distribution data: using CTEQ distribution function. 
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FIG. 5: Gluon distribution as a function of rapidity at RHIC energy. 

energy distribution as a function of rapidity at RHIC energy. We can see clearly majority of 
partons move forward or backward in number and energy, which are mostly valence quarks. 

On the other hand, fig. [7] shows the number of gluons produced and the energy 
distribution as a function of rapidity at LHC energy. The average energy per gluon is 
^f/^y" ~ 5-7 GeV at central rapidity region. This average energy is more than 10 times 
smaller compared to those of jets. 



RHIC: CTEQ: Q_0=1 .7GeV/c: Gluon 




rapidity: y 
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RHIC: CTEQ: Q_0=1 .7GeV/c: u-quark 
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FIG. 6: u-quark distribution as a function of rapidity at RHIC. 



LHC: CTEQ: Q_0=3.5GeV/c: Gluon 




Rapidity: y 

FIG. 7: Gluon distribution as a function of rapidity at LHC energy. 

Fig. [H] shows the number of u-quarks produced and the energy distribution as a function 
of rapidity at LHC energy. 

Fig. [9] shows the number density dN/dy at central rapidity region using the KNV 
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LHC: CTEQ: Q_0=3.5GeV/c: u-quark 




Rapidity: y 

FIG. 8: u-quark distribution as a function of rapidity at LHC. 

function as a function of the saturation momentum, Q s , where the saturation momentum 
Q s is different from the minimum momentum transfer Qq. Assuming that the distribution 
is flat inbetween -2.5 and +2.5 and falls off to zero at ±5 linearly, the number of parton 
produced can be obtained by integration over transverse momentum and rapidity. This 
number however depends on the saturation momentum, Q s , which is introduced in CGC 
scheme and estimated to be 1 GeV/c < Q s < 2 GeV/c at RHIC energy. The total number 
of initial partons is about 5400 for the Q s = 1.7 GeV/c and 4200 for the Q s = 1.5 GeV/c. 
We can see that the gluon density at central rapidity region at Qo = 1.7 GeV/c is compared 
to that of KNV results with Q s = 1.7 GeV/c at RHIC energy and that at Q = 3.5 GeV/c 
compared to Q s = 2.3 GeV/c at LHC energy. 

Fig. [10] and [11] show the rapidity distribution of quarks and antiquarks at RHIC and at 
LHC, respectively. We have more d-quarks than u-quarks since the Au or Pb have more 
neutrons than protons. The difference between a quark and its antiquark is the valence 
quark. 

Fig. [12] and [13] give the energy and pr distributions of produced gluons and u-quarks at 
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FIG. 9: Number distribution of KNV distribution. 




FIG. 10: Number distribution of quarks and antiquarks. 



LHC. 



Fig. [TH and [15] give the energy and pr distributions of produced gluons and u-quarks at 
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FIG. 11: Rapidity distribution of quarks and antiquaries at LHC. 




FIG. 12: Energy and pr distribution of gluons at LHC. 



RHIC. 



14 



LHC: CTEQ; Q_0=3.5GeV/c: u-quark 




GeV 



FIG. 13: Energy and pr distribution of quarks and antiquarks at LHC. 



RHIC: CTEQ: Q_0=1 JGeV/c: Gluon 




GeV 



FIG. 14: Energy and px distribution of gluons at RHIC. 
IV. SUMMARY AND CONCLUSIONS 

We showed total energy liberated from heavy ion collisions at RHIC and LHC energy 
and estimated the minimum momentum transfer, Q . We calculated the number of par- 
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RHIC: CTEQ: Q_0=1 .7GeV/c: u-quark 




GeV 



FIG. 15: Energy and pt distribution of quarks and antiquarks at RHIC. 



tons in detail produced from the collisions. We presented the rapidity distribution, energy 
distribution and pt distribution. 

These distributions can be used as initial states to evolve the system. The partonic 



transport theory will be ideal to study the evolution of the system 14 



15 



22 



23 



24 



25| as 



a function of time. We further note that these initial distributions could provide the initial 
condition for the hydrodynamic equation of motions. 
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